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Abstract 

Magnesium diboride MgB2 has been an extraordinarily "hot" research topic in the two years since 
the discovery of superconductivity below Tc « 40K in this compound in 2001. A large part of the 
excitement is due to the vast amount of unusual properties originating from the involvement in super- 
conductivity of two sets of bands that are of a very different nature. Because of the different dimen- 
sionality of the two sets of bands, this "two band superconductivity" leads two a complex behavior of 
the anisotropic superconducting state properties as a function of magnetic field and temperature. The 
interplay between two band superconductivity and the anisotropic mixed state properties and phase 
diagram is reviewed, with a strong focus on experimental results obtained by torque magnetometry on 
MgB2 single crystals. The different dimensionality of the two sets of bands is manifested in the upper 
critical field Hc2 anisotropy 7^^, which was found to strongly decrease with increasing temperature. 
While the angular dependence of Hc2 follows roughly the predictions of anisotropic Ginzburg-Landau 
theory AGLT, small, but systematic deviations were observed near T^. This, and the temperature 
dependent 7//, witness a breakdown of the AGLT description of MgB2, even close to T^. Theoretical 
calculations are in qualitative agreement with the observed ^h{T) dependence, but suggest a difference 
between 7// and the low field penetration depth anisotropy 7a with much lower values of 7a at low 
temperatures. Measurements of reversible torque vs angle curves in the mixed state on the one hand 
are incompatible with expected "single anisotropy behavior" . On the other hand, they indicate that 
the penetration depth anisotropy 7a in intermediate fields has to be higher than currently available 
theoretical estimates. The field dependent effective anisotropics observed in torque as well as many 
other experiments can be accounted for qualitatively by the faster depression of superconductivity in 
the more isotropic bands by the applied field. The vortex matter phase diagram is drawn, based on 
measurements of the reversible and irreversible torque in the mixed state. In the irreversible torque, a 
peak effect (PE) was observed in fields of about 0.85 iJc2- History effects of the critical current density 
in the PE region suggest that the PE signifies an order-disorder transition of vortex matter, analogous 
to findings in low as well as high superconductors. 
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1 Magnesium Diboride: an unconventional "conventional supercon- 
ductor" 



Magnesium diboride MgB2 is a binary inter-metallic compound ^ with a simple, layered structure: 
Boron atoms form a honeycomb lattice with hexagonal magnesium layers sandwiched between the B 
layers (see Fig. The material has been known since the 1950's T, but little attention was paid to 
this compound. This changed dramatically in January 2001, when it was reported that MgB2 becomes 
superconducting below Tc ~ 39K ,2^. Particularly the high transition temperature caused a lot of interest 
in this compound and it's physical properties. Just ten months after the discovery of superconductivity 
in MgB2, a comprehensive review jSl appeared, covering the first months of research about this "new" 
superconductor. The review counted 263 studies that between January and July either appeared in 
journals or were available to the community as e-prints. 

The initial excitement was driven by the high transition temperature, thought to be close to the 
limit of what can be expected for a material exhibiting superconductivity with a conventional, phonon 
mediated pairing mechanism 4,. The high of MgB2 therefore immediately raised speculations about an 
unconventional pairing mechanism. However, studies of the boron isotope effect on 5 and, e.g., the ^^B 
nuclear spin-lattice relaxation rate [Hj soon indicated a BCS type s-wave phonon-mediated mechanism 
of superconductivity. The rather conventional nature of the superconducting pairing in MgB2 is also 
evidenced by the absence of a boron isotope effect on the penetration depth ff] , in contrast to observations 
on cuprate superconductors 8 . 

On the other, MgB2 is far from being just an "all ordinary" conventional superconductor well de- 
scribed by BCS theory, as many other experiments showed. Point contact spectroscopy (Hj indicated, for 
example, an energy gap to ratio much lower than predicted by BCS theory, the results of specific heat 
measurements ^3^3 were also at odds with BCS predictions, and the total (Mg and B) isotope effect 
on Tc was found to be substantially lower than expected as well ^2] • 

The huge difference between the effects of boron and magnesium isotope substitution ^2j provided 
an early clue for a highly selective electron-phonon coupling (EPC), and the specific heat measurements 
|10|lllj even indicated the existence of a second superconducting gap. Calculations of the band-structure 
and electron-phonon interaction led to the description of MgB2 as a peculiar "two band" superconductor 
^l^^El, whose features show up in almost every conceivable measurement [THj . 

The chapter is organized as follows. In Sec. [3 the concept and origin of two band superconductivity 
in MgB2 are sketched, as well as its experimental confirmation. One of the features heavily influenced by 
the two band nature of superconductivity is the angular and temperature dependence of the upper critical 
field Hc2-, related to the coherence length. In Sec.|21 torque measurements of Hc2(T, 0) are presented and 
discussed. Section 13 also contains a brief discussion of anisotropic Ginzburg-Landau theory and basic 
experimental procedures (see Angst |17| for a more comprehensive discussion). The question of how 
exactly two band superconductivity shapes the anisotropy of the basic length scales below the (bulk) 
Hc2-, is still not completely solved. In Sec. 01 we review torque (and selected other) measurements 
in the mixed state and compare to recent theoretical concepts and predictions. Apart from two band 
superconductivity, MgB2 is also interesting in that thermal fluctuations of vortices in this superconductor 
are of a strength intermediate between the ones of the cuprate and of conventional superconductors. In 
Sec. El we present and discuss the vortex matter phase diagram, drawn based on torque measurements. 
Finally, in Sec. El we summarize our discussion and present the main conclusions. 
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Figure 1: Crystal structure of magnesium diboride MgB2. The structure consists of alternating layers of 
magnesium (large spheres) and boron (small spheres) atoms. Figure from Ref. jJH]) © copyright Nature 
Publishing group. 

2 Two band superconductivity: a text book example 

A two band scenario leading to two superconducting gaps was already proposed just two years after 
formulation of the BCS theory Two different bands in which superconductivity can exist are the 
most essential ingredient. Without any scattering of electrons by phonons from one band to the other 
one, there would be two transition temperatures. In the case of weak, but finite, interband phonon 
scattering, the lower disappears and the temperature dependence of the lower gap becomes strongly 
non-BCS. With very large phonon interband scattering, there should be still two gaps with temperature 
dependences now roughly similar to BCS predictions. The initial two-gap proposal was made in view 
of applications to s — d metals. However, while some hints of "two-gap effects" were experimentally 
thought to be detected in form of small deviations from BCS in transitions metals jJOj, no strong effects, 
unequivocally attributable to two superconducting gaps, were found. It was pointed out |U that, in 
order to observe considerable two-band effects, it was not enough to have two overlapping bands crossing 
the Fermi level. An additional requirement was that the bands in question have very different physical 
origins, such as covalent vs metallic-type bonding, which is more likely to occur in compounds. 

The most convincing historical demonstration of two-gap superconductivity was made by tunnelling 
experiments in Nb doped SrTiOs 22i, where clear double peaks in the tunnelling conductance were found. 
Two- or multi- band models were also proposed for the cuprate HTcSC (see, e.g., Kresin et al. [231)) but 
strong electron-electron correlations and magnetism make those materials quite difficult to tract with 
first principles calculations (24j, and unequivocal experimental evidence for any observable effects being 
due to multi band effects is scarce. It was only after the discovery of superconductivity in magnesium 
diboride that strong two-band effects were seen by a large number of different experimental techniques. 
What is it that really makes MgB2 so special? 

Answering this question was helped tremendously by the fact that the relatively simple crystal struc- 
ture (see Fig.^) and the lack of complications arising from spin-dependent interactions or strong electron- 
electron correlations made accurate ah initio calculations of the electronic structure |25p26p27j and of the 
lattice dynamics |28[l29| inn] possible. In the remainder of the section we briefly sketch this answer, about 
which by now a general consensus developed. 

Magnesium atoms in MgB2 are effectively ionized i-e., they serve mainly as electron donors, 
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Figure 2: Fermi surface of MgB2, consisting of four separate sheets: the tubular networks in the center 
of the figure are the two vr sheets, the two P centered (approximate) cyhnders the a sheets (see text). 
The coloring/shading indicates the wave vector dependent energy gap of MgB2 without any impurity 
scattering at 4K, with the coding and the distribution of gap values given besides the main figure. On 
the nearly-cylindrical a sheets, the gap value is ~ 6.5 — 7.5meV, while on the 3D vr sheets it is ~ 1 — 3meV. 
Figure from Ref. ^^1, © copyright Nature Publishing group. 

similar to electron doping of graphite by intercalation. Like carbon atoms in graphite, boron atoms are 
sp^ hybridized, with three of the four [HJI valence electrons tied up in strong covalent a {px,y) bonds 
lying in plane, while the fourth electron is in nonbonding vr (pz) states, which are delocalized (like the vr 
bonds in benzene). 

Apart from doping the boron sublattice with electrons, Mg atoms, respectively Mg^+ ions, have 
another crucial impact on the electronic structure of MgB2 : Since the Mg^"*" attractive potential is felt 
more strongly by the electrons in the vr states than by the ones in the a states, magnesium induces charge 
transfer from the a to the vr states, accompanied by the corresponding relative shift of the a and vr band 
energies. 

This leads the a bands of MgB2 to be incompletely filled, in contrast to graphite. The covalent 
a bonds are "driven metallic", i.e. the corresponding carriers become mobile, by Mg induced charge 
transfer. The calculated Fermi surface (PS), shown in Fig. [21 consists of four separate sheets. Two 
sheets, derived from the boron vr (pz) bands, are (electronlike and holelike) tubular networks. Since the 
smaller B— B in plane distance is compensated by a smaller {ppir vs ppcr) hopping jjjj, they are virtually 
isotropic and clearly three-dimensional (3D). The other two sheets, derived from boron a {px,y) bands are 
holelike (approximate) cylinders with axis P— A. The a sheets show very little dispersion along P— A and 
are thus nearly two-dimensional (2D), which could be expected from the covalent directed nature of the 
associated intraplanar B— B bonds. Correspondingly, a high anisotropy of the Fermi velocity, averaged 
over the a sheets, was calculated [2^]: {v^^fj) a / {v1) a — 46, where denotes the average over the a 

Fermi sheets. The quasi two-dimensionality of the a bands also leads to a sizeable contribution to the 
density of states (DOS) at the Fermi level — despite the very small hole doping level |26j . In fact, the 
contributions of the 2D cr bands to the total DOS at the PS is almost equal in size to the one of the 3D vr 
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bands (44% vs 56% TF ) . This feature of two sets of bands with different dimensionahty, but comparable 
contribution to the DOS at the Fermi level, is very peculiar. 

For phonon mediated superconductivity, the phonon "band structure" and the EPC are just as im- 
portant as the details of the structure of the electronic subsystem. The layered structure of MgB2 is 
reflected in the phonon dispersion curves as well, by weak dispersion of the optical branches along F — A, 
and by anisotropic acoustic branches. Of the four distinct branches with non zero energy at the zone 
center F, three were found to have no sizeable coupling to the electronic subsystem. This is not the 
case, however, for the doubly degenerate mode with E2g symmetry at F, which involves only in-plane 
(hexagon distorting) movement of boron atoms and is strongly anharmonic j3UII32j : This movement leads 
to large variations in the overlap of orbitals of neighboring B atoms, obviously modulating bond ener- 
gies, particularly of the covalent a bonds. Detailed zone-center frozen phonon calculations indeed 
find a coupling of the E2g phonon mode to the a bands, which is not only very large, but also strongly 
non-linear, suggesting even pairing via two-phonon exchange |13| . This anomalously large coupling is 
evidenced, for example, by an extremely large linewidth of the E2g mode [33|. It is important to note 
that the large and nonlinear coupling, the anharmonicity, as well as the small carrier density in the a 
bands are intimately related and have competing effects on Tc |35j . 

With such an anisotropic (sheet— or, more general, /c— dependent) EPC, the superconducting pairing 
(and order parameter, gap) should be anisotropic as well, at least in the clean limit. An analysis of pairing, 
decomposed into the four Fermi sheet contributions, led Liu et al. ^3] to propose the two-gap scenario 
for MgB2 : Similar anisotropic Fermi velocities and EPC of the two a Fermi sheets, and of the two vr 
Fermi sheets, respectively, allow the reduction to an effective two-band model. Allowing for different 
order parameters for each of these two {a and vr) effective bands increases the effective coupling constant 
relevant for superconductivity and is, together with the extremely strong coupling between one phonon 
mode and the a bands, responsible for the high of MgB2 jSH]- Later numerical calculations using fully 
anisotropic Eliashberg theory |15pi4j qualitatively confirmed the simpler effective two band model by Liu 
et al. . The gap distribution of gap values over the Fermi surface resulting from these calculations is shown 
in Fig. 121 Note that these results apply for hypothetical MgB2 without any impurity scattering ^35, . In 
the 2D a band, the superconducting gap-to- Tc ratio is slightly enhanced with respect to BCS theory, 
while the gap in the 3D vr band is about two to three times smaller, and the corresponding gap-to- Tc 
ratio is substantially below BCS predictions. 

In order to have two superconducting gaps staying open up to the bulk Tc, it is not only necessary to 
have a sizeable coupling between the two sets of bands by scattering of carriers by phonons. In addition, 
scattering of carriers between the two sets of bands by impurities needs to be very low; otherwise, any 
gap differences would be "washed out" by impurity scattering, accompanied by a Tc lowering In 
MgB2, the a — TT impurity scattering is fortunately very small even in rather dirty samples because a 
and vr bands are formed from different local orbitals orthogonal on an atomic scale and therefore cr — vr 
hybridization is very small [SHI- other hand, the intraband variation of the gaps as calculated for 

the clean case and shown in Fig. [21 was argued to be completely washed out in any realistic sample 
due to intraband impurity scattering [H5] . 

The two band two gap scenario is meanwhile well verified experimentally, from a large number of 
measurements both on polycrystalline material and single crystals: The band structure calculations are 
in agreement with de Haas-van Alphen j39ll4Uj and angle resolved photoemission spectroscopy re- 
sults, and the strong coupling between the E2g phonon mode and the a bands is evidenced by inelastic 
neutron scattering [30j and Raman measurements: Spectroscopic measurements with various tech- 
niques, including point contact spectroscopy [l ^ HH l Ulj. scanning tunnelling spectroscopy |lHllinill7| llH] . 
photoemission spectroscopy jl^, break junction ("S— I— S") tunnelling spectroscopy [HOj, and Raman 
scattering |511l52j observe directly two gaps. Moreover, the two different gaps and their location in k 
space on the a and vr Fermi sheets have been directly observed recently by angle-resolved photoemission 
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spectroscopy (ARPES) [HSl- The unusual temperature dependence of (zero field) thermodynamic proper- 
ties of the superconducting state, as seen by bulk probes, particularly specific heat pT7]lllll54j and thermal 
conductivity [SSI, ^^^o be accounted for quantitatively ^l|51J|5n] with two-gap superconductivity, 
as can the incomplete isotope effect jSlElEI- Due to the existence of two, coupled, superconducting 
condensates with different energy and length scales as well as very different anisotropics, a variety of 
interesting and unusual effects can be expected in the vortex state in non zero magnetic field, as will be 
discussed in the next sections, which focus on results by torque magnetometry. 

In summary, magnesium diboride is peculiar in that two sets of bands with very different average gap 
values contribute to the DOS at the FS with similar strengths, and that they have different "dimension- 
alities" (or more accurately very different Fermi velocity anisotropics) and different coupling strengths to 
the phonon subsystem. For the actual observation of two band superconductivity, it is also crucial that 
there is some coupling, but at the same time extremely low impurity scattering between the two sets of 
bands. Two band superconductivity in MgB2 leads to very peculiar effects showing up in a vast range of 
experiments. Experimental results in zero field can be well explained quantitatively within the two gap 
model. 



3 Upper critical field Hc2- breakdown of anisotropic Ginzburg-Landau 
theory 

In MgB2, as we have seen in Sec. |21 the crystal structure as well as the electronic and phononic band 
structure are all far from isotropic. This should lead to anisotropic superconducting state properties as 
well, for example to an anisotropic upper critical field Hc2 ■ Early measurements on polycrystalline or thin 
film MgB2 j57ll58ll59ll6UII61ll^ as well as electrical transport measurements on single crystals |62ll63ll64ll65j 
indeed found such anisotropy, but with a very wide span of values reported. 

The upper critical field is the maximum field for which superconductivity in the mixed state can 
persist. The mixed state becomes energetically unfavorable when the vortex cores overlap too strongly, 
since the order parameter is suppressed in the cores. Considering the isotropic case, the vortex core 
radius is given by the coherence length ^, the typical length scale for variations of the order parameter, 
which (ignoring the difference between GL and BCS coherence lengths) at zero temperature is related to 
the superconducting gap by 

Co = /i7;FAA(0), (1) 

where vp is an average velocity at the Fermi surface. The upper critical field is given by 

H,2 = '^o/2<2 (2) 

with the flux quantum $o = 2.07 x 10"'^ Gcm^ [66] . 

Anisotropic superconductors are usually treated within the phenomenological anisotropic Ginzburg- 
Landau theory (AGLT), which is obtained from the isotropic Ginzburg-Landau (GL) theory through 
replacing the effective mass 37^ m* in the GL free energy functional |681l66j by an effective mass tensor, 
with values m*, ml, and m* along the principal axes [HHlEnil^- In the most usual case of uniaxial 
anisotropy, all anisotropy is incorporated into the anisotropy parameter 7 = {m*/m'^f^y/'^, a simple 
constant. From the free energy variation it follows that 

7 = {m*Jm:,)'/' = XJXab = U/ic = h!^'/h!^ = 7^- (3) 
Here, A is the penetration depth, the typical length scale for variations of the (local) magnetic fleld, given 

by 

A- = 1^^, (4) 
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where Ug is the superfluid density (density of superconducting carriers), m* is the (band averaged) 
effective mass tensor [HZj, e is the electron charge and c the speed of hght. 

When the field H is not applied parallel or perpendicular to the c-axis or the a6-plane, scaling relations 
apply, for example the AGLT angular dependence of the upper critical field is given by 

H^^Ho) = h!^/{cos' e + sin2 ehlfi^ ^ Hl^/e{e), (5) 

where 6 is the angle between the field and the c-axis of the sample |72| l73j . Since effective masses [HZI 
are temperature and field independent, at least as long as polaronic effects and the thermal expansion of 
the lattice can be neglected, the anisotropy parameter 7 is assumed to be constant as well. 

In materials with complex Fermi surfaces, the situation may not be as simple as outlined above. 
For example, slight temperature dependences of the Hc2 anisotropy, which we call 7// hereafter, were 
observed in some materials, such as NbSe2 |74j . and explained by microscopic theories taking into ac- 
count non-local effects [TniIZS]- These theories show that the temperature dependence of 7 cannot be 
attributed to an anisotropy of the band effective mass tensor, unless it is also wave- vector dependent. An 
anisotropic energy gap, caused by an anisotropic attractive electron-electron interaction, however, can 
lead to variations of 7 with T. 

In MgB2, the situation is particularly complex, due to the two sets of bands: the tt band has a nearly 
isotropic Fermi velocity and a small superconducting gap, which would lead to an isotropic coherence 
length of the order of ~ 50 nm. The highly anisotropic a band with a larger gap would have a 
coherence length anisotropy of about 6, with ^o''^'' ~ 10 nm. Which set of bands controls the bulk upper 
critical field and it's anisotropy 7jy? 

The torque magnetometry results reviewed in this section (JZHZHUZnUHOl, as well as several other 
experimental l81"82"8:-{"."84"85"8 e)llS7U88j and theoretical 89'¥0 , 91797,"^ studies indicate that the answer 
strongly depends on temperature: the T dependence of 7// in MgB2 is much more pronounced than ever 
observed before in any compound. 

In order to precisely determine the angle dependence and the anisotropy of the upper critical field of 
MgB2, we carried out torque measurements in fields up to 90kOe on several MgB2 single crystals with 
volumes between about 4 x 10~^ mm^ and 2 x 10~^ mm^ Crystals of MgB2 were grown with a high 
pressure cubic anvil technique, described in detail by Karpinski et al. j78 |l95j . In brief, a mixture of Mg 
and B was put into a BN container and a pressure of 30— 35kbar was applied. Growth runs consisted of 
heating during Ih up to the maximum temperature of 1700— 1800°C, keeping the temperature for 1— 3h 
and then cooling to room temperature during 1—2 h. Flat crystals were up to 0.8 x 0.6 x 0.05 mm'^ in size, 
with sharp transitions at about 38— 39K. The crystals selected for torque measurements had — 38 K: 
See Fig. 1 in Ref. [HHI for a low field magnetization curve in the region of the superconducting transition 
of one of the crystals used for the torque study (crystal B 

The torque r = mxi?~mxff, where m is the magnetic moment of the sample, was recorded as a 
function of the angle 9 between applied field and the c-axis of the crystal for various fixed temperatures 
and fields or as a function of the magnetic field at fixed angles and temperatures. Measurements were 
performed on miniaturized piezoresistive cantilevers specifically designed for torque magnetometry |97j . 
For measurements close to in fields up to 14kOe, a non-commercial magnetometer with very high 
sensitivity was used. For the measurements in the mixed state described in Sec. |1J a vortex-shaking 
process by an additional small oscillating magnetic field perpendicular to the main field was employed to 
speed up the relaxation of the vortex lattice into it's reversible state PSIIU^ . The crystal labelled A was 
measured in this system. The other two crystals were measured in a wider range of temperatures down 
to 15 K in a Quantum Design PPMS with torque option and a maximum field of 90kOe. 

Four examples of torque vs angle curves are given in Fig. |31 Panels a) and b) correspond to measure- 
ments at 22 K. For fields nearly parallel to the c-axis, both curves are flat, apart from a small background 
visible in panel b). Only when H is nearly parallel to the 06— plane there is an appreciable torque signal. 
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Figure 3: Torque r vs. angle data of an MgB2 single crystal (B), under four different conditions. The 
raw data have been antisymmetrized around 90deg in order to subtract a symmetric background. 6c2 
indicates the angle for which the applied field is the upper critical field. The schematic drawing in c) 
shows the definition of the angle 0. After Ref. j79j . 



The curve can be interpreted in a straight-forward way: for H parallel to the c— axis the sample is in the 
normal state, while for H parallel to the a6— plane it is in the superconducting state. The crossover angle 
&c2 between the normal and the superconducting state is the angle for which the fixed applied field is the 
upper critical field. From the existence of both superconducting and normal angular regions follows im- 
mediately that i7j2(22K) < 19kOe and 85kOe < i/j2^(22K). In panel c), on the other hand, the crystal 
is seen to be in the superconducting state for all values of the angle 0, and therefore 4kOe < lfj2(34K). 
Finally, the data in panel d) show only a small background contribution [100) . Therefore, the crystal is 
here in the normal state for any 9, and we have hI^^{UK) < 14kOe. 

From figure 01 we therefore have two limitations for the upper critical field anisotropy 'jh, without 
any detailed Hc2 criterion, and without any model fits : 

jh{22K) > ^ ~ 4.47; 7i^(34K) < ^ = 3.5. (6) 

With the limitations of Eq. Q, the upper critical field anisotropy of MgB2 cannot he temperature 
independent] As a consequence, the standard anisotropic Ginzhurg- Landau theory AGLT does not apply 
for MgB2. The deviation is rather pronounced, within a change of temperature of about 0.3 Tc, 
changes, at least, by a fifth of it's value. Concerning the real values of and Hc2{0), we only have an 
estimation so far and a more detailed analysis is necessary. 
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Although, in the hght of the discussion of Fig. |21 it is clear that AGLT, with it's effective mass 
anisotropy model, is not able the describe the data measured at different temperatures consistently, the 
detailed analysis of the angle dependence of Hc2 is based on AGLT. We will show that as long as we stay 
at a fixed temperature well below Tc, AGLT is able to describe Hc2{0) remarkably well jlOlj . Although the 
location of 6c2-, for example in Fig. [3^), seems clear at first sight, this clarity disappears, when examining 
the transition region in a scale necessary for the precise determination of 9c2 (see Fig. 1 in Ref. [ZZI ) . For 
a strict analysis, it is necessary to take into account that the transition at Hc2 is rounded by fluctuations 
of the order parameter around the value minimizing the free energy. 

To assess the importance of thermal fluctuations in a superconductor, the so-called Ginzburg number 
Gi = ^{'ykBTc/H^{0)^^fj{0))'^ |lfl2j can be evaluated. While this dimensionless quantity is of the order of 
10"^° to 10-"^ in low Tc superconductors [Mllinil, it can become larger than 10 ^ in HTcSC cuprates 
|l()5j . due to the higher Tc, the very short coherence length, and the pronounced anisotropy. In the case 
of MgB2, with parameters we will obtain from Hc2{T,6), we can estimate Gi ~ 10~^: The importance 
of fluctuations in MgB2 is halfway in between the one of fluctuations in HTcSC cuprates and the one of 
fluctuations in conventional low Tc compounds. 

Fluctuations are difficult to tract theoretically in full generality. In zero and low fields, fluctuations 
in high Tc superconductors have been successfully described with the "3DXY model" of critical fluctu- 
ations (see Schneider and Singer jlU6j . and references therein) In sufficiently high magnetic fields, on 
the other hand, a different approach, the so-called "lowest Landau level" (LLL) approximation, was 
used successfully, both in high and low Tc compounds, to describe the effects of fluctuations around 
Hc2 [66 llin7|rin8[lin9|llinj . The basic physical idea 107 is that in a uniform field H\\c, the fluctuating 
Cooper pairs are moving in quantized Landau orbits characterized by kz and n. The mean field transi- 
tion temperature Tc2{H) is the temperature, at which the n = (i.e., the lowest) Landau level becomes 
stable, since Cooper pairs do not exist in higher fields, in the mean field approximation. Close to (and 
above) Tc2{H), the lowest Landau level is expected to dominate the fluctuation contributions. For a 
field high enough, H > i^LLL, and temperatures close to Tc2{H), the approximation is well fulfilled, but 
upon approaching, along Hc2(T), the zero field critical temperature Tc, the contributions from Landau 
levels with n > can no longer be neglected (see Lawrie |llll for a theoretical discussion of the limits of 
the LLL approximation). At the same time, the critical region is approached, and for low enough fields 
critical fluctuations dominate. 

In the case of MgB2 jll2j . we expect the LLL approximation to be valid, except possibly very near 
Tc, and the fluctuations to be three dimensional in nature jl7j . The fluctuation magnetization M of a 
3D system in the vicinity of the transition temperature Tc2{H) in the LLL approximation is given by a 
universal function F of T — Tc2{H) |109j : 

M ^ ( A{T-Tc2) \ 

H i7i/3 \^ (rF)2/3 y ' 

where A is a material constant related to the Ginzburg number. Eq. ((7| applies for an isotropic su- 
perconductor as well as for an anisotropic one. However, in the last case the angular dependence of 
F(A(T — Tc2)/{TH)'^/^) should be taken into account. The equation can also be applied to the torque 
of an anisotropic superconductor, by using a scaling transformation, as done by Buzdin and Feinberg for 
arbitrary isotropic scaling laws |113j . We find that the rescaled torque signal 

P = -Te^/'^{e) /[sm9coseH^/^{l - l/jjj)T^/'^'^ (8) 

with e{9) = (cos^ 9 + sin^ ^ llu)^^"^ i iii the LLL approximation, a universal function of T — Tc2 with a 
fixed value P(0) at T = Tc2{fl) [TTIj . 
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Figure 4: Selected rescaled torque P [see Eq. ^] vs angle 9 curves measured on crystal B at 22 K. From 
left to right, curves shown were measured in H = 24, 30, 40, 55, 70, 80, and 85kOe. Two criteria used 
for the determination of 9c2 are indicated by horizontal dotted lines (see text). The 9c2 values obtained 
employing the lower criterion are shown in Fig. [S] below. After Ref . [f^ • 



Taking into account the F(0) value for the theoretical dependence of the universal function for a 
3D system jllOj . we can estimate that for a volume of the sample of 8 x 10~^mm^ the rescaled torque 
P reaches at T = Tc2{H) a value of about 2 x 10"^'^ dyncmOe~^'''^K~^/'^. Figure |1] shows a selection 
of torque vs angle curves rescaled according to Eq. Q, in different fields at 22 K. The crossing of 
the P(^) dependence for each field with the line of the constant value of 2 x 10"^'' dyncmOe~^'''^K~^/^ 
determines the Hc2{9) dependence as it is shown in Fig. El below. Since there are some uncertainties 
in the determination of the volume of a small crystal, we used for comparison an additional criterion of 
6 X 10"^'' dyncmOe~^''^K~^/^, i.e., three times higher. The two criteria used are indicated by dashed lines 
in Fig. 01 It is important to stress that the results obtained depend not very sensitively on the criterion 
chosen and it will be shown later (see Fig. that with both criteria defined here, we get very similar 
temperature dependences of Hc2 and 7. Additional t[H) measurements at fixed angle give Hc2{9) values 
corresponding well to those from t{9) measurements. 

Since, for the purpose of the determination of Hc2{9), the transformation of the scaling variable is 
not necessary, no full scaling analysis was performed. From the resulting Hc2{9) curve, the anisotropy 
parameter is then extracted by an analysis with AGLT, i.e., Eq. @. We note that in the rescaling of 
the torque according to Eq. ©, the target parameter 7jy is used, which is obtained only later with Eq. ©. 
Therefore, scaling analysis and determination of with Eq. © had to be performed iteratively in order 
to self consistently find 7//. However, the 9c2{H) and Hc2{9) points obtained with the scaling analysis 
depend not very strongly on the value of used in Eq. (jH}, and the iterative procedure converges rather 
fast. 

FigureEl where all 9c2{H) and Hc2{9) points, obtained from the scaling analysis of the data measured 
at 22 K on crystal B, are plotted, shows that the data obtained employing the LLL scaling analysis are 
very well described by the AGLT formula [Eq. ©], with an anisotropy parameter = 5.1. Alternative 
curves from fits with Eq. Q, but using 7// = 4 or 6, clearly fail to describe the data. The figure 
also shows data points obtained from crystal C |115j . The irreversible properties of crystals B and C 
are rather different, i.e., they have a different defect structure. The good agreement both in value and 
angular dependence of Hc2 of crystals B and C that is observable in Fig. [21 indicates that such differences 
in the defect structure do not influence the upper critical field much, at least in the region between 22 
and 34 K jll5j . and therefore cannot influence our conclusion of a T dependent Hc2 anisotropy. We found 
the angular dependence of Hc2 to be well described by Eq. © at all temperatures not too close to Tc, 
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Figure 5: Upper critical field i?c2 vs angle 9 at 22 K. The data points were extracted with a scaling 
analysis from t{0) data measured in fixed field (o from crystal B, ■ from crystal C, cf. Fig. 0} and from 
t{H) data measured at fixed angle (c>), employing the lower of the two criteria described in the text. 
Also shown are theoretical curves according to AGLT, with anisotropy parameters = 5.1 (full curve), 
4 (dashed curve), and 6 (dotted curve). After Refs. [771175] . 

but with a temperature dependent anisotropy parameter jh{T). 

The first theoretical Hc2{T,6) calculation based on a two band model (using the clean limit, i.e., 
ignoring both inter- and intraband scattering) 90 confirmed the temperature dependent 7//(T), but 
indicated systematic deviations of the angular Hc2 dependence more pronounced at low temperatures 
jll(ij . predicting that the AGLT angular dependence should be followed well near Tc. This is at variance 
with the experimental results, as we will discuss now. 

The upper critical field Hc2i0), determined from t{9) curves measured in various fields at 33 K ~ 
0.87 Tc, is shown in Fig. |HK). By definition, r is for H\\c and H\\ab, and small for field directions 
close. This is why there are no data close to 0° and 90°. The best fit of Eq. to the data is indicated 
by the full line. Small, but systematic deviations can be seen, especially when plotting the difference 
between experimental data and best fit vs 6 (inset): at 0.87 Tc, Hc2{0) is not (accurately) described by 
Eq. ©. Deviations from Eq. © were not observed at lower T (cf. Fig. [SJ- Of course, in principle our 
experimental limitation of fields up to 90 kOe may prevent the observation of deviations from Eq. Q at 
low temperatures (< 20 K). On the other hand, an independent investigation jSHl of Hc2{0), determined 
from electrical transport measurements jll7j in fields up to 280 kOe, found no deviations of the angular 
Hc2 dependence from the prediction of Eq. Q, even at 5K. 

Deviations most pronounced in the region of 0.9-0.95 Tc were also found in a recent calculation j93j 
assuming zero interband, but high intraband impurity scattering {dirty limit, see also Gurevich |91j). 
In order to compare the deviations observed in torque experiments to the predictions of Golubov and 
Koshelev we calculated "AGLT deviations" a{e) = Hc2{e)/H^^^{9). For H^^ = 4.75 kOe and 
•jH = 3.47, a{9) has form and magnitude [Fig. Et)] very similar to deviation functions for calculated [93j 
Hc2{9) (full line) at the same temperature ,118^ . Although the theoretically predicted 93_ — 4.86 is 
higher than our data indicate, the similarity of the AGLT deviation suggests that (intraband) impurity 
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Figure 6: a) Upper critical field Hc2 vs angle of a MgB2 single crystal (A), at 0.87 Tc (symbols). 
Free fit of AGLT angular dependence H^^^{e) [Eq. ©] (full line; = 4.56 kOe, 7^ = 3.28) shows 
clear systematic deviations, highlighted in the inset. AGLT dependence with the same parameters as 
used in panel b) is also shown (dotted line, -ffjg = 4.75 kOe, = 3.47). b) "AGLT deviation" a{6) = 
Hc2{0)/H^^^{9) (see text) of the data of panel a) vs angle 9. Full line: calculation of Golubov and 
Koshelev for the same reduced temperature. After Ref. [80] . 



scattering cannot be neglected in theoretical descriptions of Hc2 in MgB2 , even in single crystals of highest 
quality. 

The upper critical fields parallel and perpendicular to the layers obtained with the scaling analysis 
and Eq. © are shown in Fig. C^). Results obtained from two crystals in two magnetometers and using 
three different Hc2 criteria are depicted as different symbols. It can be seen that the values of i^Jg (^) 
and H^}.2^{T) do not depend very sensitively on the exact Hc2 criterion used. The T dependence of H^}.2 is 
in agreement with (isotropic) calculations by Helfand and Werthamer jll9j . with H^^2{^) — 31kOe. On 
the other hand, H^}.2^{T) exhibits a slight positive curvature near Tc (due to the lack of low temperature 
data, for H^}.2^{<S) only an estimation 180 kOe < H^}.2^{Q) < 230 kOe can be given). The positive curvature 
of //Jf (T) may also account for the positive curvature of Hc2 observed in earlier measurements of 

polycrystalline material Z - Correspondingly, the anisotropy 7^ = H^}.2^ /H^}.2 systematically decreases 
with increasing temperature [Fig. IZt>)]) as we already expected from Fig. 01 A change of the criterion 
used for the determination of 6c2 leads to small shifts of the magnitude of 7//, but the temperature 
dependence is always the same. The highest upper critical field anisotropy "fu — 6 was obtained at 15 K, 
the lowest anisotropy — 2.8 at 35 K. From Fig. [7| we estimate 7//(Tc) = 2.3 — 2.7, while at zero 
temperature, may become as large as 8. 

Since the time our torque measurements were performed, many other authors also observed a tem- 
perature dependence oi^u in MgB2 as well 83,84,81,82,85,86,87,88,.12i],- Before making any quanti- 
tative comparisons, we note that electrical transport measurements j83ll84j yield too high values of -ffj2 ' 
due to surface superconductivity ("i^cs") |HZj or more complicated effects jHJ- All bulk measurements 
(torque [77|[T^ . magnetization P ^I85|l87; i88j . thermal conductivity [HI], and specific heat |HH1IH7] ) agree 
well on the H^}.!^^^) dependence and value. Concerning /fj2^(T), and consequently jh{T), however, re- 
ported values differ from each other. It is difficult to say exactly what part of the discrepancies is to be 
attributed to the difference in experimental methods and procedures, and what part originates from the 
different samples used. The deviations of the angular dependence of Hc2 from Eq. © near Tc lead to a 
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Figure 7: a) Upper critical field Hc2 vs. temperature T. Open symbols correspond to H \\ ab, full 
symbols to H \\ c, from fits of Hc2{0) data to Eq. ((S)). Up triangles are from measurements on sample 
A (with 6c2 determined with a simple "straight line crossing" criterion) and squares (circles) are from 
measurements on sample B, using two criteria of constant rescaled torque P as defined in Eq. 0. Full 
lines are theoretical curves calculated by Miranovic et al. IHOj (see also text). The dashed line indicates 
"i?^", a crossover field above which superconductivity in the vr bands is drastically suppressed, b) 
Temperature dependence of the upper critical field anisotropy 'yn = -f^j2^/-f^i2 ' determined from fits of 
-ffc2(^) to Eq. (jSJ. Symbol forms correspond to the same criteria/crystals as in panel a). The full line is 
again from the theoretical calculation by Miranovic et al. jHUl- After Ref. j79j . 



slight underestimation of /-^c2 from our torque measurements in this region, most pronounced for 
T ~ 0.9 — 0.95 Tc. This may explain at least partly the distribution over a wider temperature range of the 
change in anisotropy compared to experiments measuring directly H^^2^{T) and //Jgl^)- This problem 
should not affect the data at lower temperature, where the values of H^}.2^{T) and 7//(T) obtained in this 
work is closer to the upper end of the spread of experimental reports. We note that the limited field 
range (< 90kOe) in the experiment inevitably reduced the accuracy of our analysis in the case of the 
lowest temperature data shown in Fig. d but since no systematic deviations from Eq. @ are indicated 
in our and other measurements, we think the points shown are still reliable. It is possible that a small 
number of stacking faults may affect the H^}.2' value in this region, as discussed by Angst ;17|. However, 
a comparison of Hc2{0,T) of crystals from the same source, but grown with slightly varied procedures 
and pronouncedly different pinning properties yielded no significant differences (cf. Fig. [5] and Angst et 
al. [Zni)- On the other hand, the low temperature upper critical field of crystals from different sources 
can be more different, which may be due to significantly different intraband impurity scattering rates in 
the a and the vr band [lOj. 

The anisotropic upper critical field of MgB2 has been calculated by several authors using different 
approaches [5^11101 l^inHES]- Apart from the report by Posazhennikova et al. [55], which was based 
on an anisotropic s wave model, all calculations were based on the two band model (see Sec. EJ- The 
calculations agree qualitatively with each other as well as with experiments. For comparison, the curves 
calculated by Miranovic at al. jHUl are also shown in Fig. d Concerning the details, however, there are 
variations between the theoretical reports, as in the case of the experimental ones. For example, using 
similar input values from band structure and electron phonon coupling calculations (cf. SeclJ}, values of 
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7iT(0) of 5.7 inni, 4.6 [22] and about 6.1 ^ were obtained. 

A major difference between the calculations of Refs. |9()|l92j and the ones of Refs. j9 H I9,S j is that the 
former neglect all impurity scattering, while the latter use an approach |121j based on the a and the vr 
band being in the (intraband) dirty limit. It was demonstrated that in this case the upper critical fields, 
their anisotropy, and the temperature dependence of these quantities depend sensitively on the intraband 
scattering rates in the vr and in the a band |91j . Indications that this may indeed the case in MgB2 
are given for example by a much lower anisotropy with a decreased temperature dependence in carbon 
substituted MgB2 |122j and also by the wide spread of values of Hc2 in various samples in early reports jSj . 
Furthermore, the differences seen even in the case of high quality single crystals [lUI and the agreement 
of experimental results with the calculations by Golubov and Koshelev [HH] concerning the deviations of 
the angular dependence of Hc2 from Eq. © (see above) strongly suggest that intraband scattering is 
important in all present samples of MgB2. We note that there remains the problem that experiments, 
e.g. de Haas- van Alphen [HHlHOj . magneto-optics |12,Sj . and scanning tunnelling spectroscopy |1H] indicate 
the vr bands to be in the (moderately) dirty limit, but the a bands to be rather in the (moderately) clean 
limit whereas the comparison to experiments of the calculations carried out in the dirty limit indicates 
higher scattering rates in the a bands |124j . 

Physically, it seems clear from a comparison of Fig. d and the electronic anisotropy in the a bands 
(~ 6.8^, see SecHJ that the low temperature upper critical field and its anisotropy is dominated by these 
bands with the stronger superconductivity (larger gap). In the hypothetical situation of no coupling 
between a and tt bands, there would be two different upper critical fields for these bands with H^2 nearly 
isotropic and much smaller than H^2i due to the smaller superconducting gap in the vr bands [cf. Eqs. 

and ((21)] . Due to the coupling by interband pairing and quasiparticle exchange, H^2 becomes a mere 
crossover field |125j with superconductivity in the vr bands induced by the a bands superconductivity [3H1 
even above this field. This crossover field is indicated in Fig. [3 by a dashed line. This issue is discussed 
in more detail in Sec. ID The vr bands retain some influence on the bulk upper critical field due to 
coupling, but at low temperature this influence is not very large for the small interband pairing strength 
characteristical for MgB2 |92j . 

In the limit T — > Tc, on the other hand, thermal excitations lead to a thorough mixing of states over the 
whole Fermi surface, and the microscopic derivation of the anisotropic Ginzburg-Landau equations j71j 
is expected to hold, even if this is not the case at lower temperatures (cf. Ref. [El)- This leads to an 
anisotropy parameter being a simple average over the Fermi surface of the electronic anisotropy as a 
function of the wave vector, either weighted by the gap value (clean case, ~ 2.6 in MgB2) or evenly 
weighted (dirty case, ~ 1.2 in MgB2) [7Tlll26ll9L)j . The pronounced temperature dependence of the upper 
critical field anisotropy of MgB2 is, then, a crossover between the above two cases. 

Although anisotropic Ginzburg-Landau theory should hold for T — > Tc, Fig. [3 clearly indicates that 
the variation of jh with temperature is the strongest close to Tc. Furthermore, also the deviations of 
the angular dependence of Hc2 from AGLT predictions is observed to be particularly strong near Tc 
[see Fig. Therefore, in MgB2, AGLT seems to have an exceptionally limited range of applicability 
indeed. According to the calculations by Golubov and Koshelev |9,'-{j . the AGLT description is indeed 
limited to temperatures less than 2% away from Tc. Since such a small temperature region is difficult to 
accurately probe experimentally, one may indeed say that there is a complete breakdown of anisotropic 
Ginzburg-Landau theory in MgB2. 

In this small AGLT region very close to Tc, Eq. ^ should hold, i.e. with the definition 7^ = Xc/^ab, 
we should have ^h{Tc) = 7a (^c) — 2.6. At lower temperatures, where AGLT no longer holds, however, 
there is no reason why ^\{T) should be equal to 7h(T). Indeed, calculations of ")x{T) {in the limit of very 
low fields) 126 127 found the penetration depth anisotropy to increase with increasing temperature, from 
close to 1 at OK to the "AGLT anisotropy" at Tc, exactly opposite to the behavior o1^h{T). Furthermore, 
Golubov |127j predicted a strong deviation of the angular scaling of the penetration depth from the AGLT 
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dependence [Eq. ©], essentially due to the very small a bands c axis plasma frequency |128j . This has 
not yet been experimentally confirmed, but as we have seen above, such deviations have been found in the 
case of the upper critical field. The so far only report on direct measurements of the low field penetration 
depth on single crystals jl29j uses a method that is not very sensitive to the absolute anisotropy value 
at zero temperature, but concludes that there is very little temperature dependence of 7^- This clearly 
implies 7^ 7^ Jh for most temperatures. Such a difference of anisotropics should influence the mixed 
state properties, particularly the torque far below Hc2- 

Indeed, the field delimiting the mixed state from below, the lower critical field Hd, is mainly controlled 
by the penetration depth. From the energy needed to create an isolated vortex, it follows (within 

anisotropic Ginzburg-Landau theory) H^}^ = $o/47rA^j x ln{Xab/Cab)- Early measurements of the lower 
critical field of MgB2, performed mostly on polycrystalline material, yielded values between 250 and 
480 Oe. Meanwhile a few measurements on single crystals were reported. A direct measurement on 
single crystal by SQUID magnetometry gave values of ffJ^(O) ~ 380 Oe and H^}i^{0) ~ 270 Oe, with an 
anisotropy of 1.4 |2H1- A later direct measurement with a vibrating sample magnetometer yielded the 
surprisingly high values of -ffJi (0) ~ 2800 Oe, i/J"''(0) ~ 1350 Oe, with an anisotropy of 2.1 jl30j . Finally, 
direct measurements of i?J(5K) ~ 700 Oe (0) ~ 220 Oe were reported |85j . It should be noted that 
direct measurements of the lower critical field of single crystals are very difficult, because (i) only the 
penetration field Hp, where the first vortex enters the sample, is obtained, which can be higher than Hd 
for example due to geometrical barriers |131j , (ii) demagnetizing effects lead to an underestimation of Hp 
if not taken into account properly, particularly in thin crystals for H\\c. From our torque measurements, 
we obtain at low temperatures for different samples values in the region of Hp ~ 350 — 500 Oe, with 
little dependence on the field direction. This is compatible with the early measurements and would 
be compatible with a low penetration depth anisotropy [12611127] . However, in light of the difficulties 
associated with the direct determination of Hd and the wide spread of values reported for single crystals, 
we have to conclude that the situation regarding the lower critical field of MgB2 is much less clear than the 
situation regarding the upper critical field. For example there is no report about the angular dependence 
of Hd, where deviations from the AGLT description may occur like in the case of Hc2 discussed above. 
Often the values of Hd, as well as those of Hc2 are deduced indirectly from the magnetization in the 
mixed state, for example on the basis of the London model (see, e.g., Zehetmayer et al. j85|)- The mixed 
state Hd < H < Hc2 is discussed in the next section. Sec. |3 

4 Anisotropies for H < Hc2- there are still open questions, but "the 
smaller the field the more isotropic" 

To obtain information about the anisotropy (and A and in the mixed state Hd ^ H ^ Hc2 (existing in 
superconductors with high enough k = X/S,), an analysis based on the "London model" 66, was often used. 
This method was for example used with great success in the case of the high Tc cuprate superconductors 
(see, e.g., Refs. |132lll33ll34ll98ll35j ). where the direct measurement of the Hc2 anisotropy is very difficult 
due to the blurring out of the transition by fluctuations (cf. Sec. E}- 

The basic assumption of London model is not high k and H ^ Hc2 , but a constant magnitude of the 
order parameter throughout most of the space. Only in this case the variation of supercurrents and local 
magnetization is "simple" in most of the region and the small rest of the region (essentially the cores) 
can be treated in semi-phenomenological cut-off. This gets translated to "no core overlap", restricting 
the London regime to fields sufficiently below Hc2- In the London limit, Hd <^ H <^ Hc2 and k ^ 1, 
calculations by Kogan and coworkers jl36lll37lll38j of the magnetization components longitudinal and 
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Figure 8: Torque r of MgB2 single crystal A vs angle at 32 K in 5 kOe. a) Raw data before shaking 
(open circles). The full line is a fit of Eq. Q to the average of data collected with angle 9 increasing 
and decreasing. For 6 = 90deg, a peak in the irreversible torque can be seen (see text), b) Raw data 
collected with a "shaking" procedure (full squares). The full line is a fit of Eq. Q to the "shaked data". 
After Refs. [TTllTTn. 



perpendicular to the field lead to a torque given by |134j 



V'^oH 1 ^^ sin 261, / r]H^^ 



IM 7^ h (9) 



where e{6) = (cos^ 6* + siv? 6 /^I^Y/"^ , V is the sample volume, and 7efT is an effective anisotropy, equal 
to the upper critical field anisotropy if Eq. © is valid. We note that Eq. @ is only approximately 
valid, since the London approach ignores the contribution of the vortex cores. From studying a more 
elaborate variational model it was, however, found that Eq. Q still is a good approximation in a range 
of intermediate fields, although with renormalized parameter ry and renormalized penetration depth |139| 

Eni. 

On the other hand, Eq. Q describes the reversible equilibrium torque, whereas in low fields and at 
low temperatures, the torque usually contains an irreversible part, that depends, on the direction of the 
field change, or in constant field, on the direction of the change of the angle between field and c axis of the 
sample. The irreversible part is caused by pinning due to unavoidable point-like defects, due to extended 
defects such as stacking faults, and in highly anisotropic materials even due to the layered structure itself 
(see Sec. El for a discussion of irreversible properties). 

The combination of these effects significantly restricts the region where Eq. Q describes measured 
torque data. In fact, it is quite possible to have a material where there is no region where Eq. Q is 
strictly applicable. In the past, the equilibrium torque was often approximated by just arithmetically 
averaging the measured torque for H/6 increasing and decreasing. However, such a procedure was 
shown ,98, to lead to systematic errors in the deduction of the anisotropy. A better way to obtain the 
true equilibrium torque is to speed up the relaxation of the vortex lattice into it's reversible state by an 
additional small oscillating magnetic field perpendicular to the main field ("vortex shaking") [^H]. Note 
that this "shaking" procedure is different from performing minor hysteresis loops (cf. Sec.|SJ. A detailed 
theoretical explanation of this speed up of the relaxation to the equilibrium by the shaking procedure 
was recently given by Brandt and Mikitik jl42j . 

Figure El shows an example of one of the t{6) curves measured on crystal A for an analysis with Eq. 
©. The comparison of panels a) and b) demonstrates the importance of the vortex shaking technique for 
an accurate analysis. Whereas the curves obtained without shaking [panel a)] show a substantial degree 
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Figure 9: Field dependence of the effective anisotropy 7eff (open symbols, see text), obtained from an 
analysis of reversible torque t{6) data measured on crystal A with Eq. ©. After Ref. [77]. Also shown 
are two anisotropics determined (see Fig. E3) on crystal D at 8K (t) and 37 K (a). 



of irreversibility, this irreversibility is removed after shaking [panel b)] and the data are well described 
by Eq. Q, with an anisotropy parameter 7eff = 3.2. The "unshaked" data of panel a), on the other 
hand, can be described by Eq. © as well, but with a slightly higher least-squares deviation, and with 
a higher anisotropy parameter, 7eff = 3.4. Although the difference in 7eff is not very large and by the 
eye the fitting curves in the two panels are difficult to distinguish, the data seem to indicate that in 
MgB2 , the presence of irreversibility leads to an overestimation of the anisotropy |141j , similar to earlier 
findings in the case of YBa2Cu307_5 j98j . In Fig. a pronounced peak is visible in the irreversible 
torque for = 90deg, reminiscent of the well known "intrinsic pinning" of highly anisotropic cuprate 
superconductors. This feature is discussed by Angst and coauthors [751IT7] . 

A number of additional torque vs angle curves were measured on crystal A not far away from in 
relatively low fields. Checking first the applicability of Eq. Q, with respect to the limitations sketched 
above, we note first that all data were measured well above Hd and that the curves were well reversible 
after shaking jl43j . We believe that due to the low Ginzburg number of MgB2 (as compared to the 
cuprates), critical fiuctuations should not appreciably influence the torque in the region of temperatures 
measured. The remaining limitation is the influence of Hc2 fluctuations and vortex core contributions 
in high flelds. The upper fleld limitation of the applicability of Eq. @ is not any well deflned line in 
the H-T diagram, but a comparison with the results of Hao and Clem |14flj indicates that for the data 
shown below (Fig. EJ, Eq. © should be a good approximation, with the possible exceptions of the points 
H = lOkOe/T = 32K and = 5kOe/r = 34 K. 

The result of the analysis of reversible torque t{9) data, measured at on crystal A at temperatures 
between 27 and 36 K in fields from 1 to lOkOe, with Eq. Q |77j is presented in Fig. |2l (open symbols). 
There is a number of points worth noting: 

• The effective anisotropy 7efr is field dependent, increasing nearly linearly from 7efr ~ 2 in zero field 
to 3.7 in lOkOe, again in disagreement with the anisotropic Ginzburg Landau theory. 

• No clear temperature dependence of 7efr is visible in the range of temperatures of 27 — 36 K and 
fields below lOkOe. 

• The effective anisotropy 7cfT obtained from the analysis with Eq. Q is different from the upper 
critical field anisotropy (see Sec. 
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Before proceeding with the discussion of these results, we compare with other torque measurements 
analyzed with Eq. Q that were reported in the literature. One report of an analysis of torque data 
measured at 10 K on a small MgB2 crystal with Eq. © concluded that in that case, 7eff was field 
independent jl44j . However, much of the data used in that analysis clearly is measured in fields above 
the applicability region of Eq. Q, partly even above H^}2' analysis of the remaining data can 

support a field dependent ^eS |141j . Torque measurements performed on larger crystals were reported by 
Zehetmayer et al. |85'. In this case, data spanning the range of temperatures from 5 to 30 K and measured 
in fields of 5, 10, and 20kOe were collected. From the report it is not clear which field values were used 
at what temperatures, and what was the extent of irreversibility, but the values reported are all in the 
range between 3.5 and 4.5, monotonously decreasing with increasing temperature. Finally, Perkins et 
al. |13()j report one measurement of the transverse magnetization with a vibrating sample magnetometer. 
At 25 K and lOkOe, they obtained 7efr « 2.1. 

The different values reported are not completely inconsistent, if the effective anisotropy 7efr rises with 
H in low fields, but eventually saturates to values approaching the upper critical field anisotropy in 
fields of the order of 10 kOe. For a thorough analysis, it would be necessary to have for each measurement 
the exact information about field and temperature condition (to assess the limitation H <C Hc2) as well 
as the extent of irreversibility (given that irreversibility leads to an overestimation of 7eff, as discussed 
above). 

However, given the difference between the penetration depth and coherence length anisotropics |145j 
suggested by theoretical calculations ^126. .90) . it is instructive to check the meaning of 7efr and the 
anisotropics appearing in Eq. Q in detail. This is underlined by the difference between 7efj and the 
upper critical field anisotropy ^h- The anisotropy parameter enters Eq. Q twice via e{9). In a first 
approximation |781ll7j . the appearance outside of the logarithm (determined by the interaction between 
vortices) can be thought of as due to the penetration depth anisotropy 7a, while the appearance in the 
logarithm (involving the vortex core cutoff) is linked to the coherence length anisotropy (see also 
Ref. jnS])- A field dependence of the effective anisotropy 7eff similar as shown in Fig. IHl follows directly 
from a penetration depth anisotropy much smaller than the upper critical field (or coherence length) 
anisotropy, even if they are both field independent [TsllTTj . Recently, the torque in the London limit has 
been rederived by Kogan jl46j for the general case of 7a 7^ "fH- The generalized formula is |146j : 



y$oi^ 1 \ sin 261 



(10) 



where eA,_ff(^) = {cos^ 9 + siv? 6 ^^/'^ |147j . 7a = Xc/Xab, IH = iab/^c (•••)' denotes differentiation 
with respect to the angle 9, V is the sample volume, and e = 2.718 . . . 

Based on Eq. (fTU]) and the calculations of Kogan and Miranovic [12611^0] . Kogan predicted a sign 
reversal of the torque at low temperatures in MgB2. We measured additional torque vs angle curves at 
low temperatures and in relatively low fields to experimentally check this prediction. 

Figure ITUk l shows a t{9) curve measured (on a different crystal, crystal D) in the mixed state close 
to Tc — 38.5 K. Near Tc, the difference between 7a and is small, in agreement with theoretical 
predictions j901ll46j . The common anisotropy value 7 = 2 also is in good agreement with the data 
collected on crystal A (see Fig. El. The t{9) curve measured at low T [Fig. ^Ja)] has the same sign as the 
one measured close to i.e., there is no sign change as expected jl46j for 7a <C 7h- For 7a moderately 
lower than 7//, Kogan |146j predicted a sign change only in an angular region close to DOdeg, illustrated 
with a dashed line in Fig. llOb ). Such a partial sign change is also not observed, the maximum angular 
region where it could occur given by the irreversible region (the slight asymmetry in the irreversibility 
is due to thickness variations of the crystal). Comparing the data with curves calculated according to 
Kogan jl46j . with 11^2 = 30kOe, 7^ = 6 (see Sec. ^ and various 7a, we conclude that 7a has to be at 
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Figure 10: a) Torque r vs angle measured close to T^- Dashed line: r(0) calculated by Kogan |146j 
with 7a = 7// = 2. b) T vs 6* at low temperature. Dashed line: t(G) calculated |146j with 7^ = 2, 7^^ = 6. 
After Ref. I^. 

least 2.6, considerably higher than currently available theoretical estimates |126(I127] . Alternatively, if 
7iy in 0.2 T is much smaller than in ~ ¥1^2 I145j . the absence of a sign reversal is compatible with 
smaller 7^. However, we should mention that the best description of the data with Eq, (|10|) is given by 
7a ~ IH ~ 3.3. This value is also included in the summary of Fig. El (t). This low temperature value is 
higher than the values obtained on crystal A in the same field close to Tc, indicating that the anisotropy 
in fixed field is no longer temperature independent over a larger T interval. Concerning the experimental 
check of the predictions based on Eq. (jlOj) . we note that no total or partial sign reversal was found in r(0) 
curves measured on various crystals at low temperature in fields between 300 Oe and 15kOe. However, 
for the curves measured in the lowest fields, a small partial sign change may be hidden by irreversibiliy 
remaining after shaking. 

We think that the discrepancy between the torque data (Fig. 1111)) and the predictions by Kogan |146j 
is due to one or both of the following reasons: 

• The calculations of Kogan and Golubov jl261ll27] are valid for the limit of very low fields, whereas 
7h was calculated for the upper critical field, i.e., in very high fields (cf. Sec. EJ. Instead of 7^ and 
the coherence length anisotropy 7// being very different and field independent, they could very well 
be similar or equal, both rising with increasing field. Furthermore, for the case of arbitrary field, it 
might be more appropriate to directly start with a two band model. 

• Both Eq. Q and Eq. (|10j) are based on the London model. In a wide region of the mixed state 
in MgB2, the "no core overlap" condition of the London model seems highly doubtful due to the 
"giant vortex cores" in the vr band |48j . 

We first discuss the second possible reason. As mentioned already in Sec. |21 a hypothetical vr bands 
superconductor with no coupling to the a bands but the same superconducting gap A^r would have a 
very low upper critical field: -^^2(0) ~ 3kOe. This follows directly from the low superconducting gap 
in the vr bands, using Eqs. and (|2}. Due to the near isotropic nature of the vr bands, no significant 
anisotropy is expected for U^i- '^'^ small gap and low isotropic upper critical field corresponds a 
large (direction-independent) coherence length ^,7(0) ~ 30 nm, setting also the radius of the vortex cores. 
Hence, for a pure "vr bands superconductor", the limitation to the London model given by "no core 
overlap" or H <^ Hc2 translates (for all field directions) to ^ 3 kOe. 
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The coupling between vr and a bands by interband pairing, quasiparticle transfer etc. transforms H^2 
from an upper critical field to a mere crossover in the middle of the mixed state of MgB2. However, 
even as a broad crossover, may be expected to set a new field scale significantly affecting the mixed 
state properties. Furthermore, the whole vortex structure should be affected, with each vortex having 
two cores, one in the a bands and another one with much larger radius in the vr bands. For both points 
there is considerable experimental support. 

The two superconducting gaps are for example resolved in point contact spectroscopy measurements, 
as mentioned in Sec. I2j The feature associated with the smaller (tt) gap was shown to disappear by 
application of a magnetic field of the order of lOkOe, independent of the field direction. The H^2 
crossover is also evidenced in a abrupt change of the field dependence of the specific heat coefficient 
7(iJ) = liniT^o C{H)/T (where C is the specific heat) on single crystals jl25j . The coefficient 7 is a 
(bulk) probe of the quasiparticle density of states. As far as vortex core overlap is negligible, each vortex 
is creating the same number of quasiparticles in its core, leading to j{H) = ^^H/Hc2, where 7n is the 
normal state specific heat coefficient. The unusual ^{H) dependence in MgB2 can be well explained by 
an isotropic H^2 ~ 3 — 4kOe, above which the vortex cores in the vr band overlap so strongly that the 
contribution of the vr bands to 7 reaches its normal state value 7^ ^125j (see also Nakai et al. jl48j ) . Similar 
arguments apply for the field dependence of the thermal conductivity, which is a probe of delocalized 
quasiparticles jHJ. These arguments are strongly supported by the direct observation of large vortex 
cores in the vr bands with scanning tunnelling spectroscopy (STS) In these measurements, huge core 
radii of about 50 nm are observed, which using Eq. © would correspond to H'^2 ~ 1.3 kOe. Consequently, 
a very significant core overlap is found already in fields of 2 kOe. On the other hand, the vortex lattice 
in the vr bands was still visible in fields as high as 5 kOe, which proves that superconductivity remains in 
the vr bands due to coupling with the a bands even above H'^2- This means that the order parameter in 
the vr bands is not simply gone, but still spatially varying. The difference in the values of H'^2 as deduced 
from different experiments is simply a consequence of its nature as a crossover. 

The mixed state of a two band superconductor was recently considered theoretically by Koshelev and 
Golubov 11241. based on the same approach as the upper critical field calculations of Refs. j^lHSI and 
neglecting variations of the magnetic field. Koshelev and Golubov |124j pointed out that the behavior of 
the weak band (vr bands) partial density of states (as seen by STS and relevant for the field dependent 
specific heat and thermal conductivity) can be rather different from the behavior of the weak band order 
parameter. On the one hand, the crossover in the (spatially maximal) order parameter field dependence 
below and above the weak band "upper critical field" is much smoother than the one of the (spatially 
averaged) density of states. On the other hand, the length scale for the variation of the weak band order 
parameter is closer to the strong band {a bands) length scale than would be expected from the length 
scale of the variation of the weak band partial density of states (i.e., from the large vr band vortex cores 
observed by STS). With regard to the applicability of the London model, the order parameter is the 
more relevant quantity, since it is the order parameter that is coupled to the supercurrent density, and 
thus to the local field variation around the vortex center and to the interaction between vortices |66j . 
From the relatively smooth variation of the maximum order parameter in the weak band we conclude 
that there is still a large spatial vr bands order parameter variation up to the bulk upper critical field, 
which points against the applicability of the London model in fields larger than ~ H'^2- We mention that 
recent preliminary torque measurements also point against the London model being appropriate in MgB2 
for H > H^2- 

Before considering the other possible reason for the discrepancy between torque results and the 
predictions of Kogan |146j . it is useful to list anisotropy effects in the mixed state observed by other 
methods. The anisotropy of the thermal coefficient of the specific heat ^{H) (see above) measured for 
H\\c and H\\ah can be formally defined as the ratio of fields applied perpendicular and parallel to the c 
axis that yield the same value of 7 |125j . Normally, this gives a constant value 7cp identical to the upper 
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critical field anisotropy (this follows directly from ^{H) = j^H/ Hc2)- In MgB2, however, Bouquet et 
al. jl25j found to be field dependent: It has a value of ^ 1 for < lOkOe and then rises roughly 
linearly up to the Hc2 anisotropy when the field approaches the upper critical field. This dependence is 
reminiscent of the field dependence of 7efr discussed earlier, but is the rise of 7 starting only in elevated 
fields is an essential difference. The advantage of these specific heat experiments was the sensitivity only 
to the vortex cores, without the added complication of another length scale (the penetration depth). 
Indeed, Bouquet et al. were able to separate the contributions of a and vr bands. The contribution from 
the a band has a constant anisotropy of about 6 and the contribution from the vr band is isotropic, both 
consistent with the calculated electronic anisotropy in these bands (see Sec. The effective anisotropy 
7cj, is just a weighted average. The field dependence of 7cp, sensitive to the vortex cores only, does not 
indicate a field dependence of the coherence length, but rather that phenomenological formulas should 
take into account that there are two different coherence lengths in MgB2, and ^o-. Whenever this 
is not taken into account, the anisotropy of the best "effective coherence length" needed to describe 
experimental data has to be field dependent due to the different "intrinsic upper critical fields" of the a 
and TT bands. 

An important anisotropy effect is the distortion of the vortex lattice in fields tilted away from the c 
axis, observed for example in the strongly anisotropic cuprate superconductors by small angle neutron 
scattering (SANS) |149j : the nearest neighbors fall no longer on a circle, but on an ellipse with axis 
ratio r given by (cos^ 6 + sin^ 0/7^)^ in superconductors with a single constant anisotropy parameter. In 
superconductors with no single anisotropy, we can define a lattice distortion anisotropy 7vl: 



For 7;^ 7^ 7^ the distortion is determined by the penetration depth anisotropy in the London limit 
jl5Ulll46j : 7vL = 7a • Recently, SANS measurements were carried out on MgB2 powder samples in 5 kOe 
at 2 K by Cubitt et al. ^151 . In randomly aligned powder, each grain carries a vortex lattice of its own, 
leading to a ring of diffraction instead of resolved diffraction spots. The random orientation of the grains 
and the vortex lattice distortion together lead to a broadening of this ring of diffraction. Comparing the 
measured scattered intensity distribution to the theoretically expected distribution, Cubitt et al. deduced 
an upper limit of the vortex lattice anisotropy under the given conditions: 7vl < 1-32. In case that the 
London model is applicable in MgB2 in 5kOe at 2K, this would then correspond to the upper limit of 
the penetration depth anisotropy. In light of the discussion above, the applicability of the London model 
for the experimental conditions seems questionable. Indeed, according to the calculations of Dahm and 
Schopohl jni], the vortex lattice distortion anisotropy in fields close to the upper critical field should 
be similar to ^h- The method used by Cubitt et al. |151j furthermore depends for example on model 
assumptions of the grain alignment distribution under the infiuence of the magnetic field and on the 
validity of the angular scaling of Eq. (|llj) [see also Sec. IHl- Later, SANS was also used to determine 
the vortex lattice structure and distortion on a single crystal |152j . Measurements with H\\c indicate 
a structural transformation of the vortex lattice in fields between about 5 and 9kOe, which is likely 
associated with the suppression of superconductivity in the tt bands above H^2- Measurements with 
H tilted by 45deg away of the c axis show a vortex lattice anisotropy rising superlinearly with the 
field from 7vl — 1-3 in 2kOe to roughly 3.5 in 5kOe. Differences between the value obtained in 2kOe 
and the value obtained from scanning tunnelling spectroscopy in the same field, but for H\\ah |152j 
may indicate a deviation of the angular dependence from Eq. (|11() . Similar deviations are discussed in 
Sec. 131 A deviation from Eq. (|llj) could also explain the discrepancy to the SANS measurement on a 
powder sample jl51j . The identification of the observed 7vl with the penetration depth anisotropy is 
still questionable. However, the extrapolated low field lattice anisotropy lim//^o 7VL(-f^) ~ 1-25 should 
be equal to the low field penetration depth anisotropy, since in low enough fields neither the vortex cores 
in the a bands nor the ones in the vr bands overlap and the London model should apply. Indeed, the 




(11) 
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value of ~ 1.25 is consistent with the low field theoretical calculations |126|I127] . 

For further quantitative analyses it will be necessary to make theoretical progress beyond the simple 
(one band) London model, building a sophisticated model of the mixed state built by vortices with two 
cores (or a "double core") as well as two sorts of superconducting carriers. Steps into this direction were 
undertaken jl481ll53lll541ll2l| . but what is lacking so far is a model of the variation of the local field 
around such "double core" vortices, particularly for the case of a large overlap of the larger of the cores. 
This case is very important, since it applies for the largest region of the mixed state of MgB2 {H > 
see Fig. 

Qualitatively, the picture emerging seems to be more clear. As in the analysis of zero field specific 
heat and penetration depth data [51], it is important to consider the contributions of both sets of bands 
on any mixed state property. Each of the two sets of bands has its own set of length scales, i.e., coherence 
length ^ and penetration depth A, the latter of which is linked to the partial superfluid density [see Eq. 
(jl))]. In the case of the vr bands these length scales are isotropic, whereas in the case of the a bands, 
they are highly anisotropic, with an anisotropy of about 6 — 7, in accordance to the band structure 
calculations (see Sec. |21). These intrinsic anisotropics of the vr and the a bands probably do not depend 
on field or temperature. Any mixed state property is then determined by the relative contributions from 
both bands and if described in terms of traditional single band models has to be assigned some effective 
anisotropy, depending on the contributions of the two sets of bands. At low temperatures and fields, 
both bands contribute roughly equally to the properties, in accordance to their similar contribution to 
the density of states at the Fermi level (see Sec.|2). Since the superconductivity is weaker in the tt bands 
(smaller gap), it gets depressed faster by perturbations such as an external field or thermal excitations. 
However, due to the coupling with the a bands, it remains superconducting until the a bands become 
also nonsuperconducting: the whole MgB2 is either in the superconducting or in the normal state. Upon 
application of an external magnetic field, superconductivity is depressed much faster in the weak vr bands, 
implying a decrease of the relative vr bands contribution to any mixed state property. Consequently any 
effective anisotropy used to describe the property has to rise with increasing field. How exactly that rise 
occurs depends on the quantity in question and the definition of the effective anisotropy. This applies 
for angle dependent torque measurements as well as for specific heat and vortex lattice distortion, as 
we have seen above. Increasing the temperature instead of the field in principle has a similar effect. 
However, additionally, thermally excited quasi-particles lead to an increased mixing of the two sets of 
bands, until, extremely close to Tc, the two band character of MgB2 is no longer important for the 
mixed state properties and conventional anisotropic Ginzburg-Landau theory (AGLT) applies (cf. Sec. 
ISJ. Below this tiny region, AGLT is no longer appropriate to describe the mixed state of MgB2, and even 
phenomenological theories (Ginzburg-Landau |153j . London, . . .) should be developed and used in a two 
band form. 



5 Peak effect and vortex matter phase diagram 

In Sec. 13 we examined anisotropy effects in the mixed state, which is bounded by the upper and lower 
critical fields (see Sec. EJ. From the study of cuprate superconductors is known that the mixed state 
phase diagram contains more transition lines than the upper and lower critical fields. Identified were for 
example a melting transition between a quasi-ordered vortex lattice, called Bragg glass, and a disordered 
vortex fluid jl55j . as well as an order-disorder transition between the Bragg glass and a highly disordered, 
glassy phase |156|ll57|ll58j . The latter transition |159j has been observed also in low Tc superconductors, 
such as NbSe2 |1601ll6lj . and even in the elemental superconductor Nb jl62j . but not in ultra-pure Nb 
crystals |lU3j . This transition is generally associated with a peak in the critical current density jc and 
pronounced history effects. 

MgB2 is, particularly concerning the importance of thermal fluctuations and the value of k = A/^, 
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Figure 11: Peak effect in tlie irreversible torque of MgB2 single crystal B: a) Torque t/H vs field H 
measured at 14 K and 6 = 77.5 deg. The direction of the field change is indicated by arrows, b) t/H vs 
angle 6 in 75kOe at 15 K. After Refs. i96jl29j- 

intermediate between the high Tc cuprates and low Tc superconductors. Studying the vortex matter 
phase diagram of MgB2 may thus help in understanding the phase diagrams of various superconductors 
in a unified way. 

In single crystals of MgB2 , a quasi-ordered vortex structure has been observed in low fields by scanning 
tunnelling spectroscopy [IS] and more recently clear Bragg peaks were observed by small angle neutron 
scattering |152j . showing that at least under some conditions a Bragg glass is the stable vortex phase. 
Since, by tuning the amount of quenched random point-like disorder, the stabilization of a highly disor- 
dered phase can always be favored, an order-disorder transition in fields below Hc2 should be observable 
in MgB2 as well, at least for certain impurity concentrations. We therefore looked for signatures of this 
transition by measuring torque as a function of field and angle. 

One of the torque vs field curves measured on crystal B is shown in Fig. Illb ). For better comparison 
with magnetization curves, t/H vs H is shown. A pronounced and sharp peak effect (PE) in fields close 
to, but clearly distinct from, Hc2, can be seen. The peak effect was observed in r vs ^ measurements 
in fixed field as well, an example is shown in Fig. Illb ). However, several features of the PE are better 
manifest in t vs H measurements at fixed angle, on which we concentrate below. 

A magnification of the PE region of the t{H) curve shown in Fig. Illb ) is presented in Fig. 112^ ). 
The peak is well pronounced and very sharp (full width at half maximum 4kOe in comparison to the 
location in fields of about 78kOe). Various characteristic fields are indicated: The maximum of the peak 
for field increasing (-friax) and decreasing (-ffrnax) branch of the hysteresis loop, and the onsets of the 

T I 

peak, Hon and Hon- The separation of the two onset fields is larger than the separation of the maximum 
fields, similar to the case of the cuprate superconductors |163| . Also indicated is the irreversibility field 
-ffirr, where the two branches of the hysteresis loops meet. The peak resembles qualitatively the peaks 
observed in NbSe2 jl()4pi6()j and CeRu2 jl64j . In NbSe2, there are strong indications that the peak effect 
indeed is a signature of the order-disorder phase transition |161j . 

Figure Eb) shows a measurement on crystal C with the same external conditions as in panel a). 
In this crystal, the PE is less pronounced, but still clearly discernible. It is located in slightly higher 
fields, FAiax ^ 82.4 kOe. Crystal C was grown with the same technique as crystal B, but under slightly 
different conditions. These two crystals have pronouncedly different pinning properties, particularly for 
field directions close to parallel to the layers (see Angst JZl)- The presence of the PE in two crystals 
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Figure 12: a) Magnification of tlie peak effect region of the curve sliown in Fig. Illb ) for crystal B. Also 
shown are some of the minor hysteresis loops MHL (see text) measured, labelled A— D. b) peak effect 
region of a t/ H vs. H curve measured on crystal C, under the same conditions as were used for panel 
a) . c) Magnification of MHL B shown in panel a) . MHL started from the field increasing branch of the 
full hysteresis loop are shown as dotted lines (also indicated by thick arrows), whereas those started from 
the field decreasing branch are shown as full lines. The dashed line shows the field increasing branch of 
the (normal) full hysteresis loop. After Refs. j961l79j . 

with such pronounced differences strongly indicates that the PE, or rather its underlying mechanism, is 
an intrinsic feature of MgB2. That this is indeed the case is also supported by the recent observations 
of similar peak effects on single crystalline MgB2 from other sources by electrical transport |861l87j . ac 
susceptibility with a Hall probe [5Hlll65j . and torque magnetometry T2D]. 

The peak effect shown in Fig. El is similar to the peak effect in NbSe2, where it is attributed to the 
order-disorder transition, as we mentioned above. The order-disorder transition region is accompanied by 
history effects in the critical current density (a brief discussion of this point can be found in Ref. ^ 
feature observed experimentally for example in NbSe2 |16L)ll6lj . and very useful to identify the transition. 

To investigate possible history dependences of the critical current density, we performed several minor 
hysteresis loop (MHL) measurements in and around the peak of crystal B |96j: The field is cycled up 
and down by a small amount (1.2 kOe in steps of 200 Oe) several times, ideally until the loops retrace 
each other, indicating that the vortex system reached a stable pinned state in the given field |1601I164] . 
MHL measured, within full loops, in four different regions of the PE are indicated in Fig. \T2k) (A-D). 
MHL B is magnified in Fig. I12b l: The initial branch of the MHL started from the branch of 
the full hysteresis loop (FHL) (full line indicated by arrows) clearly is below the branch of the FHL 
(thick dashed), indicating larger hysteresis. This behavior contradicts Bean's critical state model jlf)7| . 
where the hysteresis of partial hysteresis loops can never be higher than the one of the full loop. It can 
be explained by the vortex configuration on the branch of the FHL (where the MHL was started) 
having a higher jc than the vortex configuration on the branch. Repeated cycling studies furthermore 
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Figure 13: a) Variation of the irreversible torque, scaled by H sm29, vs field, at 18 K (crystal B, after 
Ref. [Zni)- The curves were measured at angles (from left to right) 9 = 15.5, 31.5, 41.5, 51.5, 56.5, 61.5, 
66.5, 71.5, 73.5, 75.5, 77.5, 80.5 and 83.5 deg. Horizontal lines indicate zero and the criterion chosen for 
the determination of the irreversibility line |166j . For the curve measured at 71.5 deg, peak maximum 
and onset are indicated (see text), b) Height of the peaks visible in a) vs angle. The peak heights both 
from zero, and from the extrapolated background signal without peak are shown. Inset: Comparison of 
the angle dependence of the peak height at different temperatures. 

indicate the meta-stable nature of the this "high field vortex configuration", at 76.5 kOe j96lll7j . In 
MHL A, similar but less pronounced effects were observed, whereas no history effects were discerned in 
studying MHL B, C, and additional MHL far below the PE. 

The findings of the detail study of minor hysteresis loops of the peak effect in MgB2 |96|ll7j can be 

I T 

summarized as follows: Between Hon and -ffrnax, pronounced history effects occur. They can be accounted 
for by the coexistence of a metastable high-field vortex configuration with high pinning and a stable low- 

t I 

field, low pinning configuration. Above //max and below i/on, no significant history effects are observed, 
indicating that there is only one vortex configuration, which is stable. The larger hysteresis width of 
MHL started from indicates pinning in the configuration stable above i/rLax to be stronger than 
pinning in the configuration stable below Hon- This is exactly what is expected for the order-disorder 
transition |1611ll7j . Particularly, the "overshoot" of an MHL over the full hysteresis loop [Fig. 112b ')] is 
difficult to explain without a "two metastable vortex configurations" scenario. Specifically, it cannot be 
explained by relaxation effects, since this should always lead to lower hysteresis of the MHL with respect 
to full hysteresis loops. 

Figure Et) shows the irreversible part At{H) = t{H^) — t{H^) of the torque, scaled by Hsin29, 
vs field, at 18 K for various angles. The scaling was chosen to minimize the angle and field dependence 
intrinsic to the torque. Since the peak is not visible at all temperatures and angles as well as in Fig. I12| 
onsets and maxima were determined from irreversible torque curves as those shown in Fig. I13b 'l. The 
-ffon was defined as the field, where the irreversible torque starts to deviate from a straight line behavior, 
as indicated in the figure for the curve measured at 71.5 deg. The onset defined in this way is close 
to Hon as indicated in Fig. 112b .) , but we note that with the determination of onsets and maxima from 
the irreversible torque, the fine details of the differences in the field increasing and decreasing branch 
of the hysteresis loops are lost. It can be seen in Fig. 113b ) that the height of the peaks varies in a 
pronounced way with the angle 9. The peak height is shown in panel b), with and without subtraction 
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Figure 14: Angle dependence of various characteristic fields at 18 K. Shown are the upper critical field 
Hc2 ({^), the irreversibility field (x) (166] . the peak maximum field -ffmax (a), and the peak onset 
field (V). Full lines are fits of the theoretical Hc2{d) dependence [Eq. ©]. Dashed lines are guides 
for the eye. Inset: Angle dependence of reduced [divided by Hc2{0), cf. Sec.|2] characteristic fields (using 
the same symbols as in the main panel). After Ref. |96j . 

of an extrapolated linear background. It indicates that the peak effect is most pronounced in the region 
of angles between about 70 and 80 deg. This angle region is independent of temperature, as shown in the 
inset of Fig. [T3b\ 

One possible explanation for this behavior is an influence of stacking faults on the peak effect. Al- 
though the presence of the peak effect, and it's location, are not affected by stacking faults, the extent 
of hysteresis may be. The difference of how pronounced the peaks of crystals B and C are [see Figs. 1111 
and 112b ')] supports such a scenario. The location in higher fields of the peak effect in crystal C indicates 
that there is less point-like disorder present in this crystal than in crystal B, if we accept the PE as 
a signature of the order-disorder transition jl561ll57lll58| . What is difficult to explain with only one 
sort of disorder varying is the smaller ratio of the hysteresis width within to the hysteresis width below 
the PE, Tirr(-f^max)/''"irr(-f^on)5 of Crystal C, as Compared to crystal B. Both the angular PE amplitude 
maximum and the difference in the ratio Tjj-j- 

(-f^max)/7"irr(-f^on) of Crystals B and C might be attributed to 
the influence of non-collective strong pinning by stacking faults in crystal B. This point was discussed in 
more detail by Angst |17j . 

The variation of the peak onsets and maxima with angle at 18 K is shown, together with Hc2{0) 
and i?irr(0) |166j . in Fig. EI These characteristic peak fields follow the angular dependence of Hc2, as 
indicated by fits to the theoretical Hc2{0) dependence according to the anisotropic Ginzburg-Landau 
theory [Eq. (jSJ of Sec. 0], while the angular scaling of the irreversibility field is less clear. This can 
be seen also in the inset, displaying the 9 dependence of the characteristic fields, reduced by the upper 
critical field. The onset field is approximately constant at about 0.8 Hc2 and the maximum field at about 
0.85 -ffc2- HiT-j- is located at about 0.9 Hc2, but seems to get slightly lower as — > 0. The good angular 
scaling of the peak onset and maximum fields supports the conclusion that they are unaffected by stacking 
faults. 
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Figure 15: Phase diagram of MgB2 single crystal at an angle of 77.5 deg between the c-axis of the crystal 
and the applied field, constructed from the curves shown in the inset. The temperature dependence of the 
characteristic fields Hc2, -f^max; and iifon is given. They mark boundaries between the normal state and the 
various phases of vortex matter, indicated in the figure and discussed in the text. The irreversibility field 
-ffirr is also shown. After Ref. [HE]- Inset: Variation of the irreversible torque, scaled by Hs'm29, vs field 
H, at 77.5 deg. The curves were measured at temperatures (from left to right) T = 27,24,22,20, 18, 16, 
and 14 K. The horizontal line near zero marks the irreversibility field criterion |166j . 



Torque measurements are not possible for H\\c or H\\ab and SQUID measurements performed with 
H\\c and H\\ab gave no indication of the peak effect. However, ac susceptibility measurements on the same 
crystals indicate the peak effect to be present both for H\\c and H\\ab, confirming that the underlying 
mechanism is a feature for all directions of H j79pi68j . This conclusion is supported also by recent reports 
of other authors, in the case of H\\c from transport data |871l86j and ac susceptibility with a local Hall 
probe |86[ll65j . in the case of H\\ab from transport data [86j . 

The phase diagram for an angle of 77.5 deg between H and the c axis of the crystal (corresponding 
roughly to the angle where the PE is most visible), is presented in Fig. 1151 In the inset, scaled irreversible 
torque vs H curves are shown. The magnitude of the peaks is reduced quickly by increasing the temper- 
ature. The main panel shows the temperature dependence of the characteristic fields. The abrupt end 
of Hojy[T) and H^i^x{T) is due to the decreased sensitivity of the magnetometer (In fact, the peak effect 
is discernible at temperatures up to ^ 27 K in measurements of the torque vs angle, such as shown in 
Fig. II lb ) (79) ) and due to the smearing of the effective pinning landscape by thermal fluctuations JJj , 
as indicated by the drop of the irreversibility line. 

In a recent report of low frequency ac susceptibility measurements |165j . the peak effect was observed 
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for H\\c at temperatures up to about 25 K, and interpreted in terms of the order-disorder transition as 
well. Pissas et al. jl65j report the transformation of the PE into a "step-like" ac susceptibility for the 
temperature interval between 25 and 27.5 K, and interpreted this as a signature of thermal melting. In 
our case, no step-like feature in the reversible torque was observed in the continuation of the PE. It 
should be emphasized, that thermal melting so far below Hc2 would be at odds j^n] with theoretical 
expectations |in5j . 

The equilibrium order-disorder transition, which corresponds to i/max t26j) is located in fields of about 
0.85 -ffc2 in crystal B and in about 0.9i7c2 in crystal C. The peak effect observed in other crystals by 
transport was reported to be located even closer to i?c2 j87ll86j . These differences are natural for a 
disorder-induced phase transition in crystals with varying degrees of disorder. Form and location of the 
PE observed in MgB2 resembles results obtained on NbSe2 single crystals with varying degrees of disorder 
|169j . but are rather different from the order-disorder transition in cuprate superconductors jl63, Lt|- 

6 Conclusions 

We have reviewed superconducting state anisotropy effects in magnesium diboride MgB2, with the main 
focus on results obtained by torque magnetometry. MgB2 is a "two band" or "two gap" superconductor: 
the superconductivity is Fermi sheet dependent in an extreme way, with a large gap (strong supercon- 
ductivity) on the (J Fermi sheets and a roughly three times smaller gap on the vr Fermi sheets. This is 
possible because of the very different coupling to the lattice vibrations in these two sets of bands and 
because of there is very low hybridization between them, resulting in a very low impurity scattering rate 
between them. Pronounced and unusual anisotropy effects originate in this two band superconductivity 
and in the very different electronic anisotropies of the two sets of bands. 

One of the unusual anisotropy effects in the superconducting state is the very pronounced tempera- 
ture dependence of the upper critical field anisotropy 7//, falling from a value close to (with some sample 
dependence) the electronic anisotropy of the a bands (~ 6.8) at zero temperature to low values of the 
order of about 2 for temperatures approaching — 39 K. There is a good qualitative agreement between 
various experiments as well as theoretical studies on this. The origin of this temperature dependence 
of is the increased "mixing" of the two sets of bands by thermally excited quasiparticles at higher 
temperatures — at low temperatures the upper critical field is mainly determined by the stronger super- 
conducting a bands, which has the higher field scale. The interplay between the different field scales, 
anisotropies, and interband mixing also leads to the observed deviations of the angular dependence of Hc2 
from conventional behavior. No general consensus has yet emerged concerning the impact of the different 
intraband impurity scattering rate on the upper critical field. Concerning the lower critical fields, the 
available experimental data is sparse, preventing the drawing of firm conclusions. 

The anisotropy effects on the upper critical field, as well as on other quantities, implies a breakdown 
of the standard (one band) anisotropic Ginzburg-Landau theory, except in a tiny (less than 1 K) region 
below Tc. Below this region, microscopic calculations have to be used, or at least a two band formulation 
of the anisotropic Ginzburg-Landau theory involving two order parameters and two sets of length scales 
^ and vr, as well as some description of the coupling between these two components. 

In the mixed state, experimentally deduced anisotropies are different from different quantities, and 
also often depend on the model used to analyze the data. However, all experimental results [at least 
those spanning an extended field region including (at low temperatures) the range of 1 — 15kOe] find 
anisotropies that are monotonically increasing as a function of the magnetic field and that are consistent 
with reaching for H — s- Hc2 or already in lower fields. In the case of torque magnetometry, describing 
the data with an equation allowing for different anisotropies of the penetration depth and the coherence 
length yields anisotropies of these length scales that are rather similar. 
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The field dependent effective anisotropies of various mixed state quantities are qualitatively well 
described by a crossover from a low field region, where the quantities are determined by similar contri- 
butions from the isotropic tt and the anisotropic a bands, to a high field region, where they are (almost) 
determined by the a bands alone. This crossover is due to the very low "intrinsic upper critical field of 
the TT bands" H^2^ order of magnitude lower than the one (in c direction) of the a bands, because of 
the weaker superconductivity (smaller gap) in the vr bands. In decoupled vr and a superconductors, there 
would be a sharp transition at -^^2' but the coupling between the two sets of band transforms this into 
a smooth crossover, which in addition may be qualitatively different for the gaps and the quasi-particle 
excitations in the vortex cores. This may also explain the qualitatively different forms of the field depen- 
dence of the anisotropy as observed for example in torque, small angle neutron scattering, and specific 
heat experiments. 

The existence of different length scales for the two sets of bands, rather directly demonstrated for the 
coherence length (for example by scanning tunnelling spectroscopy and analysis of field dependent specific 
heat data), but not yet for the penetration depth (naturally, there is only one local field distribution for 
the whole of MgB2 ) , strongly suggests that any phenomenological model should be extended to a two band 
version before being applied to the analysis of data obtained on MgB2 — otherwise obtained "best fit" 
single band values of basic length scales will be effective values that acquire a field dependence due to the 
field dependence of the relative contributions of the two sets of bands, although the bands specific length 
scales likely are not field dependent. In the case of the London model, which was used so successfully 
for the analysis of torque, small angle neutron scattering, muon spin rotation, magnetization, and many 
other data measured e.g. on cuprate superconductors, it is rather unclear whether such an extension can 
be successful. This is so because of the large overlap of the vr bands vortex cores throughout most of 
the mixed state phase diagram, which violates a basic assumption of the London model. It is clear, that 
while the puzzling mixed state anisotropy effects may be understood on a qualitative level, there are still 
open questions concerning their more quantitative description. 

From the study of cuprate superconductors is known that the mixed state can be strongly affected by 
thermal fluctuations and inevitably present material disorder, leading to phase diagrams containing more 
transition lines than the upper and lower critical fields. The material disorder leads for example to a 
destruction of the (quasi)ordered vortex lattice in fields lower than the upper critical field. This happens 
in an order-disorder transition, which was also observed in conventional superconductors, such as NbSe2. 
This transition, manifested by a peak in the critical current density with peculiar history effects is observed 
in MgB2 as well, by torque magnetometry, as well as various other measurements. MgB2 is interesting 
in terms of these additional "vortex matter" transitions as well, because of the intermediate importance 
of thermal fluctuations, compared to cuprate and conventional superconductors. A melting transition 
driven predominantly by thermal fluctuations has not been observed clearly in MgB2 so far though. 
We note that there are possible influences of the two band superconductivity on the non-equilibrium 
properties and these additional transitions of vortex matter as well. However, a direct infiuence has not 
been indicated by any experiment so far, and is thus left for future studies. 

We would like to conclude this chapter by briefly commenting on the consequences of the results 
that were reviewed for applications of MgB2. The relatively high Tc and critical current densities (not 
affected by "weak links" across boundaries), coupled with the easy and cheap processing (as compared 
to the high Tc cuprate superconductors) has created also an intense technical interest in MgB2 [T^ . 
However, the relatively large electronic anisotropy of the a bands, which inevitable becomes the relevant 
anisotropy when the magnetic fleld is increased to values depressing the superconductivity in the isotropic 
TT bands, rather degrades the possibility of granular MgB2 material to carry high critical current densities 
jc in high magnetic flelds: correspondingly, jc has been found to be a rapidly decreasing function of the 
applied magnetic fleld T. In high quality single crystals, magnetization curves are almost reversible in 
a wide region even at low temperature [see Fig. Illh .)]. It will be important to maximize the pinning 
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properties by increasing the disordered vortex phase, which is restricted to a small region below Hc2 in 
clean material (see Fig. I15|) . through the controlled introduction of uncorrelated material defects. This 
increases the processing requirements for MgB2 based materials for high field applications. However, the 
potential of strongly affecting also the anisotropy by selectively tuning the intraband impurity scattering 
rates |911ll22j may lead to applicability enhancements that cannot yet be predicted. 
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